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INTRODUCnON 


From  the  bathtub  drain  to  the  airplane  wing  to  the  mushroom 
cloud  of  a  nuclear  explosion  the  vortex  occurs  as  a  dominating  fluid 
flow  feature.  In  fact,  the  vortex,  or  "eddy"  as  it  is  sometimes  called 
in  reference  to  turbulence,  is  probably  the  most  common  and  the  most 
important  single  fluid  flow  phenomenon;  yet  it  remains  one  of  the  most 
poorly  understood.  Vortices  play  an  important  part  in  the  transition  to 
turbulence  (Reference  1)  and  in  meteorological  storms  (Reference  2) . 

Possible  applications  of  the  vortex  are  in  heat  exchangers  (Reference  3) , 
gas  separators  and  vortex  reactors  (Reference  4) ,  (Reference  5) ,  the 
Ranque-Hilsch  effect  (Reference  6),  advanced  space  propulsion  (Reference 
7) ,  (Reference  8) ,  and  power  generation  systems  (Reference  9) ,  (Reference 
10).  Additional  applications  as  well  as  an  extensive  bibliography  of 
recent  work  on  vortices  are  given  in  reference  11. 

Many  of  the  older  theoretical  analyses  of  the  vortex  utilized  highly 
restrictive  and,  in  many  cases,  unreal  assumptions  such  as  negligible 
axial  velocity  or  no  variation  of  the  tangential  velocity  in  the  axial  direc¬ 
tion.  It  is  only  through  the  use  of  such  restrictive  assumptions  that  exact 
solutions  of  the  full  Navier-Stokes  equations  can  be  obtained  for  so  com¬ 
plex  a  flow.  Recently,  Long  (Reference  12)  obtained  a  numerical  solution 
of  the  laminar  Navier-Stokes  equations  with  boundary  layer  approximations 
which  is  not  too  restrictive.  In  this  solution,  the  axial  velocity  is  of  the 
same  order  of  magnitude  as  the  tangential  velocity  for  all  radii.  It  is  noted 
by  Long,  however,  that  experimentally  the  axial  velocity  falls  off  faster 
than  the  tangential  velocity  at  large  radii.  Also  recently,  Lewellen  (Ref¬ 
erence  13)  obtained  a  series  solution  of  the  laminar  Navier-Stokes  equations . 
Newman  (Reference  14) ,  using  boundary  layer  assumptions ,  linearized 
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and  integrated  the  Navier-Stokes  equations  for  viscous  laminar  vortical 
flow  behind  a  wing  tip.  He  then  applied  the  laminar  solution  to  a  turbu¬ 
lent  vortex  by  replacing  the  kinematic  viscosity  by  an  eddy  viscosity. 

This  assumption  of  an  eddy  viscosity  which  is  a  function  of  the  circula¬ 
tion  and  kinematic  momentum  of  the  vortex  but  not  of  the  radius  or  the 
axial  coordinate  has  been  used  in  most  analyses  of  turbulent  vortices. 
Newman  (Reference  14)  states  that  the  eddy  viscosity  varies  slightly 
in  the  radial  direction  and  that  this  variation  is  about  the  same  as  in 
the  turbulent  Jet.  He  also  found  variation  of  the  eddy  viscosity  in  the 
axial  direction  in  the  wing-tip  vortex  but  attributed  this  variation  to  the 
incompleteness  of  the  rolling  up  of  the  vortex  rather  than  to  actual  variation 
of  the  turbulence  level.  Gasparek  and  Dosanjh  (Reference  15)  obtained 
experimental  velocity  measurements  in  a  vortex  behind  an  airfoil  (NACA  0009) 
and  compared  these  with  Newman's  theory. 

In  the  present  investigation,  the  turbulent  entrainment  analysis 
of  the  axl-symmetric  free  turbulent  Jet  developed  by  Paris  (Reference  16) 
is  extended  to  describe  the  free  turbulent  vortex,  and  the  resulting  theory 
is  compared  with  that  of  Newman  (Reference  14)  and  the  experimental 
data  of  Gasparek  and  Dosanjh  (Reference  15). 

The  theory  is  then  extended  to  describe  the  confined  turbulent 
vortex  and  is  compared  with  experimental  measurements  taken  in  a  sink 
type  vortex  confined  in  a  cylinder. 


EXPERIMENTAL  APPARATUS  AND  TECHNIQUES 

Apparatus 

The  apparatus  used  to  produce  and  study  the  confined  turbulent 
vortex  Is  shown  in  Figure  1.  Figures  2  and  3  show  close-ups  of  the 
5.72  inches  Inside  diameter  plexlglas  tube  with  the  vortex  core  made 
visible  by  smoke. 

The  air  entered  tangentially  at  one  end  of  the  tube  through  a 
system  of  twelve  one-inch  high  vanes  set  tangent  to  the  tube  and  placed 
symmetrically  around  its  circumference.  The  vanes  were  covered  with 
a  glass  plate  which  permitted  visual  and  photographic  observation  into 
the  end  of  the  tube . 

A  wooden  plate  with  a  1.25  Inches  diameter  hole  in  its  center 
was  located  inside  the  tube  2.50  feet  from  the  vanes.  The  plate  fitted 
flush  with  the  Inside  of  the  tube  so  that  no  airflow  around  its  edges  was 
possible.  The  air  was  pulled  through  the  hole  in  the  plate  by  a  Joy 
Axivane  fan  connected  to  a  Sola  constant  voltage  transformer.  A  variac 
was  used  to  vary  the  power  to  the  fan  and  thus  regulate  the  flow  velocity. 

The  probe  used  in  measuring  the  velocity  profiles  and  its  traversing 
mechanism  are  shown  in  Figure  4.  The  probe  was  made  of  progressively 
smaller  diameter  stainless  steel  tubing  ending  in  a  section  0.018  inch 
in  diameter.  Before  the  measurements,  the  probe  tip  was  accurately 
adjusted  to  be  exactly  on  the  axis  of  the  probe.  The  micrometer  traversing 
mechanism  was  fitted  with  a  calibrated  angle  adjustment  device  such  that 
measurements  could  be  taken  at  any  angle  to  the  axis  of  the  tube . 

Lines  parallel  to  and  in  the  same  plane  with  the  axis  of  the  tube 
were  scribed  on  each  side  of  the  tube  to  assist  in  orientation  of  the  probe. 
A  series  of  holes  spaced  2  inches  apart  was  drilled  along  one  of  the  lines 
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to  accommodate  the  probe.  All  holes  except  the  one  in  use  were  covered 
with  tape.  The  tube  was  mounted  in  its  supports  so  that  the  probe  moved 
in  the  plane  containing  the  two  lines  and  the  tube  axis.  Final  fine  adjust¬ 
ment  of  the  probe  into  this  plane  was  made  with  the  two  adjusting  screws 
at  the  base  of  the  traversing  mechanism.  The  aluminum  block  on  the  front 
of  the  traversing  mechanism  was  cut  to  conform  exactly  to  the  outside 
diameter  of  the  tube  to  insure  that  the  probe  axis  was  perpendicular  to 
the  tube  wall. 

A  strip  of  electrically  conducting  paint  was  painted  on  the  inside  of 
the  tube  along  the  line  opposite  the  holes.  Batteries  and  a  small  bulb 
were  connected  between  the  paint  and  the  probe  so  that  a  complete  circuit 
was  formed  when  the  probe  touched  the  paint.  Thus,  it  was  possible  to 
determine  the  exact  micrometer  setting  at  which  the  probe  just  touched 
the  wall  without  damage  to  the  probe.  The  angle  adjuster  was  zeroed 
visually  using  the  line  scribed  on  the  tube.  Measurements  were  made 
with  a  Kollsman  helicopter  airspeed  indicator  calibrated  against  a  Betz 
manometer. 

Techniques 

Because  of  the  very  large  velocity  gradient  in  the  radial  direction 
(the  tangential  velocity  increasing  from  zero  at  the  centerline  of  the  tube 
to  a  maximum  value  in  about  one-fourth  inch),  it  was  impossible  to  measure 
the  static  and  total  pressure  at  the  same  time  as  it  conventional.  Therefore, 
at  each  setting  of  the  micrometer,  measurements  were  taken  at  two  degrees 
of  angle  intervals  and  at  one  degree  intervals  where  necessary.  The  maxi¬ 
mum  pressure  was  then  taken  as  the  total  pressure,  and  the  minimum  pressure, 
which  was  located  approximately  90  degrees  from  the  maximum  pressure. 
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was  taken  as  the  static  pressure.  The  angle  at  which  the  total  pressure 
occured  was  taken  as  the  angle  of  the  velocity  vector. 

The  velocity  was  then  calculated  from  the  total  and  static  pressure 
measurements,  and  its  components  in  the  tangential  and  axial  directions 
were  obtained  from  the  angle  of  the  velocity  vector.  About  45  measure¬ 
ments  were  taken  between  the  wall  and  the  centerline  of  the  tube,  at 
0.20  inch  Intervals  near  the  wall  and  at  0.020  inch  intervals  near  the 
centerline . 

The  centerline  was  located  by  extrapolating  the  tangential  velocity 
curve  to  zero.  This  accounted  for  the  small  irregularities  in  the  thickness 
of  the  paint  and  in  the  tube. 

Measurements  were  taken  through  five  holes  at  the  same  variac 
setting,  the  five  holes  being  1.5,  5.5,  13.5,  21.5,  and  26.5  inches 
from  the  plate.  Also,  measurements  through  the  13.5  inch  hole  were  taken 
at  a  higher  and  a  lower  variac  setting . 


THEORETICAL  DESCRIPTION  OF  THE  FREE  TURBULENT 
SOURCE  OR  SINK  TYPE  VORTEX 


Tangential  Velocity 

The  tangential  velocity  of  an  inviscld,  potential  flow  vortex  is 
given  by 


AT 


n 


(1) 


yielding  infinite  velocity  at  (Figure  5).  However,  in  a  real 

fluid  the  velocity  discontinuity  at  will  be  smoothed  out  by 

the  action  of  viscosity. 

It  was  suggested  by  Cornish  (Reference  17)  that  it  might  be  possible 
to  obtain  the  viscous  vortex  tangential  velocity  profile  by  considering 
the  potential  flow  vortex  velocity  profile  shown  in  Figure  5  to  represent 
two  hypothetical  adjacent  Jets  in  opposite  directions .  Following  this 
analogy,  at  /i,  ^  O  there  exists  a  velocity  discontinuity  between  the 
two  Jets  which  must  be  smoothed  out  by  mutual  entrainment. 

This  turbulent  entrainment  process  was  analyzed  by  Faris  (Reference 
16)  in  his  description  of  the  axi-symmetric  free  turbulent  Jet.  Using  a 
boundary  layer  approach,  he  found  that  the  axial  velocity  in  such  a  Jet 
can  be  given  by 
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(see  Figure  6) .  The  function  is  a  function  of  / ,  where 

(i')  is  the  width  of  the  Jet,  and  is  plotted  in  Figure  7.  This  function 
is  very  closely  approximated  by 


n 

In  direct  analogy,  the  jet  on  the  right  side  of  the 

^  ^  0  Itne  in  Figure  S  entrains  flow  on  the  left  side  of  the  A  =  line 

given  by 


r 


'  JUdt 


c5l  TT/Z 


(4) 


where  now  is  the  width  of  the  zone  of  influence  on  each  Jet  on 

the  other. 

Thus,  as  is  shown  in  Figure  8,  the  entrainment  velocity  of  each 
Jet  is  opposite  in  direction  to  the  primary  velocity  of  the  other  Jet  so 
that  the  resultant  velocity  is  obtained  by  subtracting  the  entrainment 
velocity  of  each  Jet  from  the  primary  velocity  of  the  other.  This  resul¬ 
tant  velocity  profile  is  then  taken  as  the  viscous  free  turbulent  vortex 
tangential  velocity  profile  and  is  shown  by  the  dotted  curve  in  Figure  8. 

Therefore,  the  tangential  velocity  of  the  free  turbulent  vortex  is 
given  by 


AT 


Ah 
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or, 


AT 
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At) 


(5) 


As  can  be  seen  in  Figure  8,  AVi)  is  the  radius  of  the  boundary  layer 
within  which  the  effects  of  viscosity  alter  the  potential  flow  profile . 
Outside  the  boundary  layer  -  ! )  the  effects  of  viscosity  are 
negligible,  and  the  tangential  velocity  is  given  by  the  potential  flow 
expression: 


JZ_ 

^  ”  <J‘7ryi 


(6) 


Pressure 

According  to  the  usual  boundary  layer  approximations  for  vortex 
flow  (Long,  Reference  18)  the  centrifugal  force  and  the  pressure  gradient 
force  are  the  only  important  terms  in  the  radial  Navier-Stokes  equation: 


=  F  if 

Therefore,  the  pressure  is  given  by 


(7) 


(8) 
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Substitution  of  AT  from  (5)  and  (6)  and  application  of  the  boundary 
condition  at  *•  >0  ;  ^  yields  (Appendix  A) 

f=  =  -  Csh)'[FCP)  -  s.ijj  ^  i  I ^  (9) 


The  function  Is  plotted  in  Figure  9 . 

Axial  Velocity 

In  a  real  vortex  the  axial  velocity  is  in  general  not  negligible  as 
is  often  assumed.  It  is  thought  that  the  same  tangential  velocity  and 
pressure  distribution  may  exist  with  various  modes  of  the  axial  velocity 
depending  on  the  kinematic  momentum  and  boundary  conditions .  It  will 
be  assumed  here,  subject  to  subsequent  experimental  verification,  that 
in  some  cases  the  axial  velocity  is  the  same  as  that  of  the  axi-symmetric 
free  turbulent  jet.  In  this  case  the  axial  velocity  of  the  vortex  is  given 
by  (Paris,  Reference  16) 


) 


(12) 
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Following  Long  (Reference  18)  the  kinematic  momentum  transfer  In 
the  axial  direction  is  assumed  to  be  constant;  thus, 


<X2 


Ji  j  ^ =  o  . 

*0 

Substitution  of  /^J'and  from  equations  (9),  (10),  (11),  and  (12) 
yields  (Appendix  B) 


(13) 
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(14) 


or. 


constant 


=  ±  St  6 


(15) 


The  plus  sign  Is  to  be  taken  for  source  flow  and  the  minus  for  sink  flow . 
Using  (15),  the  axial  velocity  from  (11)  and  (12)  becomes 


(16) 
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Radial  Velocity 

Using  the  axial  velocity  from  (16)  and  (17),  the  radial  velocity 
may  be  obtained  by  integrating  the  continuity  equation: 
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Thus, 
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(19) 


Substitution  of  j/J"  from  (16)  and  (17)  and  application  of  the  boundary 
condition  at  >Z  -  ^  -xc  s  ^  yields  (Appendix  D) 


(20) 
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The  function  "J" { is  given  by 


T(f)  5  ^  ^/- 


^-o 


(22) 


and  is  plotted  in  Figure  10. 
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Boundary  Laver  Radius 

It  now  remains  to  determine  K(i)  .  According  to  Newman 
(Reference  14)  and  Squire  (Reference  19)  the  kinematic  viscosity  may 
be  replaced  by  a  turbulent  eddy  viscosity  S  which  is  a  function  of 
K  and  n  ,  and  Is  constant  throughout  the  vortex  to  a  good  approx¬ 
imation.  Employing  the  usual  boundary  layer  approximation  for  vortex 
flow  (Long,  Reference  18),  the  axial  Navler-Stokes  equation  then  be¬ 
comes 


>t-cr 
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(23) 


Evaluated  ~  O  ,  (23)  becomes 


Jo 


-rm 


(24) 


Since  yu.  -  0  at  ~0,  Evaluation  of  these /fe  limits  from  equations 
(9)  and  (16)  and  substitution  in  (24)  yields  (Appendix  F) 


=  +  9.0‘/ e  . 

Integrating, 
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Again,  the  plus  sign  applies  to  source  flow  and  the  minus  to  sink  flow. 

As  stated  before,  the  eddy  viscosity  ^  is  also  a  function  of  A' 
and  n  .  This  linear  dependence  of  on  ^  is  in  agreement 

with  Long  (Reference  12)  and  Gartshore  (Reference  11). 

The  constant  of  integration  in  (26)  would  be  the  radius  of  an  orifice 
at  M  "^0  if  the  boundary  layer  approximations  were  valid  at  2  -  O  . 
However,  near  the  orifice,  the  partial  derivatives  with  respect  to  5 
cannot  be  neglected,  and  (26)  does  not  apply  exactly.  Because  of  this 
and  the  existence  of  a  boundary  layer  next  to  the  wall  of  the  orifice,  the 
constant  in  (26)  will  be  slightly  less  than  the  radius  of  the  orifice. 

It  can  be  seen  from  (26)  that  in  the  case  of  source  flow  (jet  with 
circulation)  becomes  infinite  if 


r. 


<5.//  /r 


//. 


r 


(27) 

is 

is 


and  for  values  of  greater  than  0. 665  ,  Idi 

negative.  Also,  for  sink  flow  (common  sink  vortex) 
negative  for  values  of  less  than  0.665  .  Therefore , 

j  /P 

since  only  positive  values  of  are  realistic,  equation  (26)  applies 

only  to  source  flow  with  zero  or  small  circulation  as  compared  with  the 
kinematic  momentum  and  to  sink  flow  with  strong  circulation. 


Non-Dimensional  Equations  and  Certain  Characteristic  Values  from  the 
Above  Equations 

Figures  11,  12,  and  13  show  typical  variation  of  the  , 

and  profiles  in  the  axial  direction.  To  facilitate  comparisbn  with 
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experimental  data ,  the  following  values  are  obtained  from  equations  (5) , 
(9),  (16),  and  (20). 


=o.3gi>  (ir) 

J 


=  =  o.  /^O 

J  ^  i 


—  -d3.5- 


n 

AJ-  =  /.  ^  . 

/ntiuy  <B7r/\  ^ 


(28) 


ahaV 


r  ±('0. 


(29) 


AT 

''''ahaV 


A- 


r^3£>  J2_  _ 

—  0  •  S(o  /  //^\  ~f^  ~  ^  /Z 

(7f)  >  ^ 


(30) 


-  0-^3oI^F{^)'S.iTJ  ^  -0  =  ^.3^^  y2 


(31) 


W 


=  I 


yi.  -A 

30  rs  ~  y^. 


0~ tY^) 


yz 


(32) 


(33) 

A 

-  a.  7  so  'y[' 


Equations  (30),  (31),  (32),  and  (33)  are  plotted  in  Figures  14,  15, 


16,  and  17. 
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Application  to  Wlna-Tlo  Vortex 

In  the  case  of  the  wing-tip  vortex,  the  axial  velocity  outside  the 
boundary  layer  is  not  zero  as  in  the  previous  case.  Therefore,  using 
coordinates  fixed  relative  to  the  wing  with  2  positive  behind  the 
wing,  the  axial  velocity  deficiency  rather  than  the  axial  velocity  is 
assumed  to  be  a  similar  function  of  T 


^  Ad  , 

=  lA^do  ^  ^  . 


(34) 


After  Newman  (Reference  14)  all  the  wing  wake  is  assumed  to  be 
drawn  into  the  vortex,  an  assumption  which  is  not  justified  near  the 
wing  but  becomes  increasingly  accurate  as  ^  increases.  The  profile 
drag  of  the  wing  is  then  given  by  the  deficit  of  longitudinal  momentum 
within  the  vortex: 


(35) 


Assuming  -  >W  ')«  \AJ  oo  (Newman,  Reference  14),  equation  (35) 

may  be  written 


~  ^  A.  Jji  . 


0 


(36) 
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Substitution  for  from  (34)  and  integration  (Appendix  G)  yields 

(37) 


Therefore,  (34)  becomes 


\AJao  “  >w  ^ O’  y^Ko  yr^ ^  R*’  ^  ^ 


(38) 


y^. 


^  ^  \Mac  ^  TT  -  /. 


Using  the  A.-  0  limits  from  Appendix  H,  the  axial  Navier-Stokes 
equation  (24)  becomes 

r  /  Oo 

y 


C o.  /S6  ttR  ^  J 


(39) 


=■ 


Pc 


0’  7r/^(/y> 


-H. 


Integrating, 


2  -h 


constant. 


(40) 


Thus,  the  radius  of  the  boundary  layer  is  proportional  to  the  square 
root  of  £  in  the  case  of  the  wing-tip  vortex,  rather  than  to  the  first 
power  as  in  the  source  or  sink  type  vortix.  This  conclusion  was  also 
reached  by  Gasparek  and  Dosanjh  (Refenmce  15)  and  Squire  (Reference  19). 
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There  is  no  restriction  on  the  circulation  in  the  wing-tip  vortex, 
all  values  of  H  being  possible  according  to  (40) .  It  should  be  noted 
that  Newman  (Reference  14)  assumed  small  axial  velocity  deficiency  and 
small  tangential  velocity  compared  with  the  remote  axial  velocity.  Here 
only  small  axial  velocity  deficiency  has  been  assumed. 

Experimental  data  taken  in  the  wing-tip  vortex  behind  a  NACA  0009 
airfoil  in  a  wind  tunnel  by  Gasparek  and  Dosanjh  (Reference  15)  is  com¬ 
pared  with  the  theory  of  this  paper  in  Figures  18  and  19.  Though  the 
scatter  is  understandably  great,  the  agreement  is  generally  good.  As 
was  predicted,  the  agreement  between  theory  and  experiment  becomes 
greater  as  the  distance  from  the  wing  increases. 


Application  to  Axi-Svmmetric  Free  Turbulent  Tet 

If  r  -  0  ,  equation  (26)  applies  to  the  axi-symmetric  free 

turbulent  jet: 


„  /7.  yc?  6- 


-h 


constant. 


(41) 


According  to  Schlichting  (Reference  20),  for  the  axi-symmetric 

free  turbulent  Jet  is  independent  of  the  kinematic  momentum  and  is  equal 
to  0.29.  This  value  was  confirmed  by  Faris  (Reference  16).  Therefore, 
from  (41) 


/  7.  y:?  (£ 

{IP 


o- 


(42) 
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or, 


6'  =  0,0/6^  . 


(43) 


Schlichting  (Reference  20,  p.  608)  gives 


-  O.0/&  /  i/F'  . 


(44) 


The  close  agreement  between  (43)  and  (44)  confirms  the  use  of 
the  axial  Navier-Stokes  equation  to  determine  . 


THE  CONFINED  TURBULENT  SINK  TYPE  VORTEX 
Extension  of  the  Free  Vortex  Theory 

The  theory  of  the  preceding  chapter  applies  only  to  free  vortices . 

The  Introduction  of  confining  walls  produces  pronounced  changes  In 
the  vortex  as  can  be  seen  In  Figure  20  where  experimental  measurements 
of  the  tangential  velocity  of  a  vortex  confined  In  the  tube  described  in 
Chapter  II  are  compared  with  the  free  vortex  theory. 

It  was  found  In  the  preliminary  measurements  that  when  the  plate 
with  the  hole  In  It  was  removed  from  the  tube,  the  introduction  of  smoke 
showed  no  characteristic  vortex  core,  and  the  tangential  velocity  of  the 
flow  was  directly  proportional  to  the  radius  all  the  way  out  to  the  edge 
of  the  boundary  layer  on  the  wall  of  the  tube. 

Therefore,  it  was  thought  that  the  tangential  velocity  of  the  con¬ 
fined  sink  type  vortex  produced  with  the  plate  in  the  tube  might  be  described 
by  a  superposition  of  the  free  vortex  profile  and  the  linear  profile: 

r  A  » 

^  ~  -+  M  /I  .  (45, 

Experimental  Results 

The  experimental  results  of  the  measurements  through  the  five 
holes  at  one  variac  setting  are  presented  in  Figures  21A-21E.  In 
Figures  21C-21E  the  experimental  data  is  compared  with  equation  45, 
and  the  agreement  can  be  seen  to  be  quite  good  over  most  of  the  curve. 

The  circulation  P  decreases  and  jk.  increases  as  the  distance  from 
the  plate  increases .  The  tangential  velocity  at  the  two  stations  closest 
to  the  plate  (Figures  21A  and  2  IB)  is  not  a  linear  function  of  the  radius 
near  the  centerline  as  is  the  case  at  the  other  three  stations,  and. 
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therefore,  equation  45  does  not  apply  near  the  plate. 

The  tangential  velocity  curves  all  exhibit  a  characteristic  kink 
just  beyond  the  maximum  values.  The  region  of  the  kink  is  also  the 
region  of  departure  from  equation  45.  It  is  possible  that  this  may  be 
due  to  a  perturbing  effect  of  the  probe ,  for  near  the  centerline  audible 
perturbations  occured  when  the  probe  was  in  the  static  pressure  position. 
Visual  observation  with  smoke  also  revealed  that  the  presence  of  the 
probe  caused  a  slight  displacement  of  the  position  of  the  visible  core 
of  the  vortex. 

The  axial  flow  pattern  is  very  complex  as  can  be  seen  in  Figures 
22A-22E.  Throughout  most  of  the  tube  there  is  very  little  radial  flow. 
Observations  with  smoke  indicated  that  after  the  air  enters  through  the 
vanes  it  does  not  move  radially  into  the  center  but  spirals  downward 
along  the  wall  toward  the  plate.  The  smoke  forms  a  hollow  visible  core 
as  shown  in  Figures  2  and  3.  The  radial  position  of  the  smoke  core 
apparently  corresponds  to  the  second  zero  point  in  the  axial  velocity 
curves.  The  top  of  the  smoke  core  intersects  the  glass  plate  at  almost 
a  right  angle  so  that  the  core  seems  independent  of  the  incoming  flow. 
When  the  source  of  the  smoke  was  removed,  the  core  of  smoke  remained 
in  position  and  only  after  several  minutes  was  gradually  drawn  out  through 
the  plate. 

Observations  with  a  wool  tuft  indicated  that  the  flow  through  the 
hole  in  the  plate  is  down  toward  the  fan  near  the  edges  of  the  plate. 
Nearing  the  centerline  the  direction  of  the  flow  changes  to  up  toward 
the  vanes.  It  was  not  possible  to  determine  the  direction  on  the  center- 
line  because  of  the  size  of  the  tuft. 
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Therefore,  the  following  description  of  the  axial  flow  may  be  given 
as  drawn  in  Figure  23  where  the  axial  velocity  profiles  at  each  of  the 
five  stations  and  the  tube  are  drawn  to  scale.  The  air  enters  tangentially 
at  the  vanes  and  spirals  downward  along  the  wall  toward  the  plate.  At 
the  plate,  part  of  the  flow  moves  out  through  the  hole,  and  part  turns  and 
moves  upward  along  the  outside  of  the  smoke  core  toward  the  vanes . 

This  flow  turns  at  the  vanes  and  joins  the  incoming  downward  flow  along 
the  wall. 

Inside  the  smoke  core  the  flow  is  upward  toward  the  vanes  next 
to  the  core  and  downward  toward  the  plate  near  the  axis,  being  zero  on 
the  axis.  Thus,  the  axial  flow  exhibits  two  co-axial  elongated  ring  vor¬ 
tices,  one  outside  the  core  and  one  inside.  The  inside  ring  vortex  extends 
through  the  hole  in  the  plate  and  has  its  turning  point  below  the  plate. 

Observations  with  an  aural  stethescope  Indicated  that  the  flow  was 
weakly  turbulent  inside  a  radius  of  about  0 . 85  inch  and  strongly  turbulent 
inside  a  radius  of  0.30  inch,  the  transition  from  weak  to  strong  being 
gradual. 

The  results  of  the  measurements  at  different  variac  settings  are 
presented  in  Figures  24  and  25.  As  can  be  seen  the  flow  pattern  is  the 
same  at  all  three  settings,  only  the  magnitudes  of  the  velocities  being 
affected. 


CONCLUSIONS 


The  use  of  the  turbulent  wake  function  and  the  boundary 

layer  turbulent  entrainment  analysis,  which  correctly  describe  the  axi- 
symmetrlc  free  turbulent  jet,  has  been  extended  to  describe  the  free 
turbulent  source  or  sink  type  vortex.  The  agreement  between  theory 
and  experiment  in  the  wing-tip  vortex  is  good,  and  the  theory  correctly 
predicts  the  value  of  the  turbulent  eddy  viscosity  in  the  axi-symmetric 
free  turbulent  jet.  Therefore,  it  may  be  said  that  the  theory  does 
correctly  describe  the  free  turbulent  vortex. 

The  confined  turbulent  sink  type  vortex  may  be  described  by  a 
superposition  of  the  tangential  velocity  of  the  free  vortex  and  a  tangential 
velocity  directly  proportional  to  the  radius  of  the  vortex.  The  agreement 
between  this  modified  theory  and  experiment  in  a  vortex  confined  in  a 
tube  is  good  over  most  of  the  curve,  and  the  slight  departure  from  the 
theory  may  be  due  to  pertubations  caused  by  the  measuring  probe. 

Experimental  measurements  revealed  that  the  axial  velocity  in  a 
confined  vortex  is  quite  complex,  consisting  of  two  co-axial  elongated 
ring  vortices  separated  by  a  zero  velocity  region. 


APPENDICES 
Table  of  Contents 


Appendix 


Page 

Evaluation  of 

f  =f-^ck  +/rzj - 

—  26 

Evaluation  of 

-+  1^3^  A  =  - 

'o 

- 28 

Evaluation  of 

I  =  fo-^ASrf)J(w) - 

—  31 

Evalutlon  of 

—  32 

Evaluation  of 

_ 

(2? 

• —  36 

Evaluation  of 

IrX , 

,  c  ^  X 

G.  Evaluation  of  J 


37 

39 


H.  Evaluation  of 

ri  IT).  ^  ^(§fX 

for  Wing-Tip  Vortex -  4q 


24 


Appendix  A  Evaluation  of  -  f  yi.  ^  * 


Substituting  for  AT  from  (S) , 


f i(t)  - i(^) .  w^i. 


(A-1) 


Graphical  Integration  yields 


f  =  .  #  -  /, 


(A-2) 


where 


'’/k 

=  / 


=J-7£t  , 


(A-3) 


and  is  plotted  in  Figure  9 . 

Outside  the  boundary  layer,  AT  is  given  by  (6),  and  the  pressure 
by 


f  • 


+  '%(A 


/i. 


>  IT  -  I. 


(A-4) 


Integrating, 


(A-5) 
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(A-6) 


(A- 7) 


Equations  (A-2)  and  (A- 7)  must  be  equal  at  • 


and,  therefore,  since  f^O)  ^  S,  j 

=  T  -  ■ 


(A-8) 


(A-9) 


Substituting  (A-9)  in  (A-2), 


f 


e 


-  j;  If 


> 


yj. 

R 


/. 


(A- 10) 


I 


ADOendix  B  Evaluatlnn  nf  =  O. 
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From  (11)  and  (12). 


since  /cw  —  O  for  From  Appendix  C, 


/ 

fo*AjY^)Y^)  = 

^  tk 


Therefore . 

^  ‘ 

J 


=  0.0^//  {lA/"-/^^) 


Also, 


Substituting  for  p  from  (9)  and  (10) , 

^  i. 


^  J-  fs  ± 

^  Jij2T 


y  C  ^ 


(B-l) 


(B-2) 


(B-3) 


(B-4) 


(B-5) 
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J. /JL. ^  f  (b-6) 


But  the  values  of  all  three  integrals  are  constants.  Therefore,  if  is 
not  a  function  of  2  . 


And,  therefore, 


J  r‘*° 


W  =:  constant. 


But  the  kinematic  momentum  is  given  by 


or. 


/r  «  ^7r  , 


{B-9) 


(B-10) 


(B-11) 
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Therefore,  from  (B-9)  and  (B-11), 

\AJR  =  t  /7P  =  constant.  (B-12) 

The  plus  sign  is  to  be  taken  for  source  flow  and  the  minus  for  sink  flow . 


29 


I 

Appendix  C  Evaluation  of  ^  ^ 

Using  the  approximate  form  of  fr®*”  (3) ' 

I  =  TtTJ { ,  (C-1) 

Let  ^  =^/-  Tt)  .* 

=  fj}/-^)jy.  -  ^  (c-2) 


Expanding  the  cosine, 


//  '  y  ^/n. 

(~l)  TT 

y  ^ 

V 


,  . /H  Sf>t 

r-i)  TT 


^  (Hm  f  I)  ^ 


A|  =  cl 


Therefore , 

X 


/  ^  ^  I 

,»/  !  JS  ^  I 

~~~{f-y)u  ~  '^ 4 

,  i  /  ^  ^  ^  J 

^  ^yt^TTy  1^  -h  ^  ('J^J'.Ciyr^-,,)  /"  L 

/  ^  /  ^ 

”  -  y^2/UW  /*■  /  , 

or,evaluating  the  series. 


(C-3) 


(C-4) 


X  =  V-  5r/2?.  ;>ys- J  -  ,^.^3//. 


(C-5) 
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Appendix  D  Evaluation 


I  j  , 


From  (16), 


^  jr] 


(D-1) 

/ 


Therefore , 
yzM.  =  :t  cP. 


or, 


^  >u. 


+ 


(D-2) 


(D-3) 


Integrating  the  last  Integral  by  parts; 


j^2b  )  ^(“r^  ~  /4/>  ( r) 


(D-4) 
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Therefore , 


The  Integiral  is  evaluated  in  the  same  manner  as  in  Appendix  C  so 
that 

fA/f)J/'t)  =  -  rM  fS 


^-0  7r 


^19-  ^  //_ 


(D-6) 


Therefore , 


-h^^yrO-tf  -  J-  ^ 


(D-7) 


^  -  / 


where 


*•  tt 


Finally, 

^  =  "ir  Si 


(D-8) 


/X  , 


But  at  ^  ^  ,  ju.  ^  O  so  that 

aquation: 

Therefore,  since  7T<5 J  = 

//*}  =  ^ ro. m ®-'‘’> 

and  (D-8)  becomes 

(D-11) 

t  ;tr(^)*0./fTj^ 

It  is  shown  in  Appendix  E  by  application  of  L'Hospital's  rule  to 
(D-11)  that  AA.  is  indeed  zero  at  A  ~  O  . 

Outside  the  boundary  layer,  ^  =^and  the  continuity  equation 
becomes 


must  satisfy  the 


=  0 

•  (D-9) 


) 


(9, 


Therefore , 


(D-12) 
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But  (D-11)  and  (D-13)  must  be  equal  at  yz  ~  : 


t  s.  ^6  JTF'  /t  ^^3)  =• 


('2) 


since  TTO  ~  0  •  Therefore, 


('z)  =  ±  o,  ^ ^ 


(D-14) 


(D-15) 


and  (D-13)  becomes 


/. 


(D-16) 
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Appendix  E  ^valuation  of  JXyvH,  AJ.  . 

A  C) 

From  (D-U), 


J^ju.  = 

A-^o 

'  ^  ^  ^ttO- /  rr^)  -h  O.  /^?J 


(E-1) 


which  by  L'Kospital's  rule  becomes 
M^aj.  =■  t  ^ ^ 

A*9<J  Tf^o 


-  -  J/"/-  ^)‘'-3-r^^)J _ 


(E-2) 


But,  from  Appendix  D, 


,  oo  ^-f)  rr  ^  j 

T  C^)-  ”  =  -G^vCf-ir^  . 


(E-3) 


Therefore , 


A  -^o 


» 
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Appendix  F  Evaluation  of 


Iri ,  ^ 


9A  . 


From  (9) ,  since  F^O)  =  0  , 


Therefore , 


/  r  \*-  !  J/^ 

=  ^  Jb 


(F-1) 


(F-2) 


From  (16),  since  (o)  —  O  , 


j.  ffF 
/? 


Therefore , 

ff)o  ^ 


JR  \ 

/K  A  JJ J 

s,//  /r  JR 

^2  * 


(F-3) 


(F-4) 


Also,  from  (16), 


Juaj-  /  ^  *7^  •i ^ 

^  ^  ”  ~  ~p- 


(F-5) 


Similarly, 


^ 

^>2*  ”  ”  '  (F-6) 


Differentiation  of  (3)  yields 


•othat  =  -  J>r\ 


Therefore , 

(  I  ~  (  /?■»  y  .  (F-8) 


Also, 


:f 


s.  Si  /p 

/? 


(F-9) 


I 


I 

Evaluation  of 


Proceeding  as  In  Appendix  C, 


Ja 


-^/e 

*o 


C>-«.  Tr 


(G-l) 


^  -J-/0-37’'/)  =  C!.ai3.  <°-2) 


I 
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Appendix  H  Evaluation  of  ^  C  >4 


for  wing-tip  vortex. 


Form  (34),  since  ~ 

/— — £l_  \—L 


Ia4 


Therefore , 


C'-  W\  - 


o( 


/  Oo  )  /  JR 

(  TT/®  /yT'* 


Also, 


- \-a' 

o>ya.  '  K  o. /si  7r  lA/>^  J  i 


5 

a/l* 


A 


(H-1) 


(H-2) 


(H-3) 


O.  tt  z®  IA/^ 


(H-4) 
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// 

Therefore,  since  , 


which,  as  In  Appendix  F,  becomes 


The  pressure  term  in  unaltered,  belno  given  by  (F>2)  in  Appendix  F. 
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Confined  Vortex  Core,  Side  View 
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Confined  Vortex  Core,  End  View 
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Figure  22C 
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